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A randomized first-passage problem for drifted 
Brownian motion subject to hold and jump from a 

boundary 
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Abstract 

We study an inverse first-passage-time problem for Wiener process X{t) subject to 
hold and jump from a boundary c. Let be given a threshold S > A(0) > c, and a 
distribution function F on [0,-|-oo). The problem consists in finding the distribution 
of the holding time at c and the distribution of jumps from c, so that the first-passage 
time of X{t) through S has distribution F. 
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1 Introduction 

In this brief note, we generalize the inverse first-passage-time problem for Wiener process 
X with random jumps from a boundary c, studied in our paper [2], by introducing holding 
when X hits c. Diffusions with holding and jumping boundary (DHJ) were hrst introduced 
by Feller (113). and their ergodic properties were studied in [19], and in [9] in the case 
without holding. In the one dimensional case, a DHJ coincides with a diffusion in J = (c, d), 
until it hits the boundary c; when this occurs, the process waits there for an independent 
random duration following an exponential distribution with parameter /3 > 0, and then it 
makes a random jump in the interior of /, according to a given distribution z/, and starts the 
diffusion afresh. The (direct) hrst-passage time (FPT) of a one dimensional DHJ through a 
barrier S with c < X(0) < S' < d, was investigated in [IB], while in [5] the FPT was studied 
in the case without holding, i.e. (3 = J-oo (see also [20] 1. 

The drifted Brownian motion with holding and jump boundary c can be described as 
follows. Let Xn{t) = X + fit + Bt, where B is standard Brownian motion (BM), and x G 
I = [c,d]. We construct a new stochastic process X, called drifted Brownian motion with 
holding and jump boundary c: starting from x G /, at time t = 0, we set X{t) = X^{t) 
until the (random) time at which the boundary c is reached, then the process X^ is killed 
and it is continued as the new process X which holds for an independent exponential time 
at c, with the holding rate /5, then it makes a random jump from c inside I and it starts 
afresh (independently of the past history) from a point U & I which is a random variable 
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with distribution u, and then it evolves as until the boundary c is reached again, and so 
on. Indeed, in the present paper we take c = 0 and d > 0, for the sake of simplicity. 

Hitting the boundary c = 0 by the process X can be interpreted as the occurrence of a 
catastrophe with delayed effect; after holding, it consists in resetting the state of X to a new 
value, and this has various applications in the scope of diffusion models in Queueing theory 
(see e.g. m), as well as in Mathematical Finance and in Biology (see 0). As an example, 
X{t) can represent the number of individuals of a population at time t (for instance fishes 
in a little lake); whenever the population goes extinct (i.e. it reaches the level 0), after a 
waiting time, it is performed a restocking by a random amount U > 0. 

Notice that, considering X^(0) > c and c = 0, is not restrictive. For instance, let us 
consider the diffusion Z with holding and jumping boundary c > Z(0), associated to Z{t) = 
Z{0) + fit + Bt with fi > 0, Z{0) > 0, and U G (0,c). Z is similar to the Wiener-type 
neuronal model in the presence of refractoriness, considered in |H], [IE], for which the neuron 
hres when its voltage exceeds the threshold c, and after the refractoriness period, the voltage 
is reset to f/ G (0, c). We observe that the process Z can be reduced to our case; indeed, by 
using that —Bt is distributed as Hj, it follows that the distribution of the hrst hitting time 
of Z to c, when starting from Z{0) < c, is nothing but the distribution of the hrst hitting 
time of to zero, when starting from c — Z{0) >0. 

Let S' G / = (0, d) be a constant barrier, and suppose that the holding rate (3 and the 
jump distribution iy(u) = P{U < u), u G (0,5') are assigned; then, the FPT of X over S, 
when starting from 0 < a; < S', is: 

Ts{x) = rs,ji,y{x) = inf{f > 0 : X{t) > S'|X(0) = X(0) = x}. (1.1) 

The direct FPT problem for X was studied in [18], while the FPT for the process without 
holding (i.e. (3 = -1-cxd) was studied in |5]; here, we deal with the following inverse FPT 
(IFPT) problem: 

for a given barrier S E I and 0 < x < S, let F be a distribution on [0, -|-oo), 

then the IFPT problem consists in finding the holding rate (3 and the jump distribution u 

or its density g, if it exists, so that ts{x) has distribution F. 

In this case, we will say that the the pair {g, (3) is solution to the IFPT problem; notice that 
the density g has support (0, S'). 

This IFPT problem is analogous to that considered in [2] for diffusion with only jumps 
(without holding) from a boundary, and studied in [3], [1], [13] for diffusions without holding 
and without jumps, and in [T] for reflected diffusions. It has interesting applications in 
Mathematical Finance, in particular in credit risk modeling, where the FPT represents a 
default event of an obligor (see m), and in Biology (see ini)- The paper is organized 
as follows: Section 2 contains the main results; Section 3 provides some extension to more 
general diffusions and in Section 4 we report some explicit examples. 

2 Main results 

Let X be a one-dimensional diffusion on D = (0, -|-cxo) and let X the associated diffusion with 
holding and jumps from the boundary c = 0, starting from a; > 0. Following the notations 
of [IB], we denote by /5 > 0 the holding rate of X at 0, and by u the jump distribution 
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when it hits 0. Moreover, let t{x) be the first-passage time (FPT) of X through S, with 
X < S, and let hs{t,x) the density of t{x) and hs{X,x) its Laplace transform (A > 0). 
Finally, denote by G\, (A > 0) the potential (or resolvent) operator of X, which is given 
by G\i>{x) = Ip{X{t))dt, for any Borel function -0 on [0,-|-cx3), and by Gf the 

potential operator of X on iA, and killed at tdIx), i.e. the hitting time of X to 0. We recall 
the following result from [T8] . 


Theorem 2.1 Let denote the transition density of the diffusion X killed at zero; then, 
the Laplace transform of the FPT of the DHL X through S is: 


hs{\,x) 


p{X,x, S) 

Pix,s,sy 


( 2 . 1 ) 


where p is given by: 


p{X,x,y) 


p^{X,x,y) + (1 


AGfl(a;)) 


6{y) + yp^{X,u,y) 
A(1 -I- /3(z/, Gf 1)) 


( 2 . 2 ) 

□ 


Let X be BM with drift p, namely, X{t) = X^{t) = x F pit + Bt, then the transition density 
of X killed at 0 is (see [10], pg. 221): 


p^{t,x,y) 


1 

\/2'xt 


{y-tit-xf /2t 


_ /2t 


(2.3) 


By calculations analogous to those concerning the case /i = 0 (see e.g. [2T], pg. 82, [Tl|, pg. 
354, [15], pg. 288) we obtain its Laplace transform: 


P^iX,x,y) 


^y{y-x) 

v/2A + pF 


^-\y-x\y/2\+PP _ ^-\y+x\^/2\+f? 


(2.4) 


Thus, if vi^dx) = g{x)dx, by straightforward calculations and using the expression of the 
Laplace transform of the hrst hitting time of drifted BM to the barrier S (see m), (E3) 
becomes: 


X S{y) + fd jy"p^{X,x,y)g{x)dx 

p{X,x,y)=p {X,x,y) + e X J ---^ . - 

A + /3 ( 1 - /“ 


(2.5) 


where p^ is given by fl2.4p . 

Our main result is the following, where, for the sake of simplicity we limit ourselves to 
consider the IFPT problem for the DHJ X associated to BM with drift p < 0, in the special 
case when the starting point is a: = 0. 

Theorem 2.2 Let X(t) = p,t + Bt, p < 0, and let us consider the associated process X with 
holding and jumps from 0, having holding rate /9 > 0 and jump distribution ir, suppose that 
the FPT of X over S > 0 has an assigned probability density f(t) = f^(t) and denote by 
/^(A) = f^{t)dt, A > 0, the Laplace transform of f^. Then, if there exists a solution 
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(5',/9) to the IFPT problem for X, the Laplace transform 'g{X) of g, and {3 must satisfy the 
equation: 


l(\) [ (l - (a + /3 (l - g ( v'2A + + g))) 

+") (g(-V2A + g2 ) - g 

= -c) (g V2A + ) - g (V^A+T? )) 


( 2 . 6 ) 


Moreover, if p, = 0, and we require that the density g is symmetric with respect to then 


formula 


can be explicited and we get: 


or 


dW = 


/o(A) = 


/o 


AM /A2 


+ /3{l + e-^^) 


/9(l + /o(f) 


A + /3(1 + e-5^) - ^g(V^) 

Proof. Setting /^(A) = hs{X,0), by Theorem 12.11 we have 

- _ P(A,0,>^) 

^ p{X,S,S) 

By using fl2.5p and calculating the various quantities, it is easy to see that: 

/} (g(-y2A+7?)-g(y2A+7?)) 


(2.7) 


( 2 . 8 ) 


(2.9) 


P(A,0,S) = 


\/2A + 




and 


p{X, S, S) = 


X + /3 (^1 — g (^\/2X + /i^ 

1 _ g-2Sv'2A+/A2 


V2A + /i2 

P e-A^'/^+A (g(-y2AT7?)-g(y2AT7?)) 


\/2A + 




X + /3 (^1 — g ^'\/2A + /i^ 


from which (I2.6p follows. If g is symmetric with respect to I-, we have: 

?(-A) = e^^g{X). 

Then, equation fl2.7p easily follows, by taking /i = 0 in fl2.6p and by using fl2.10p . 


( 2 . 10 ) 

□ 


Remark 2.3 If there is no holding at zero (i.e. (3 = +cx)), by letting (3 —)■ +cx) in fl2.7p . we 
obtain: 

/o(t) (1 + e-"^)) 

m = \ ^ -;^gA2A -• (2.11) 


1 + /o (t) 


which coincides with equation (2.11) of [2] 
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Remark 2.4 By taking the first and second derivative with respect to A in fl2.7p and cal¬ 
culating them at A = 0, we get 

B(rs(0)) = i + 1 - 2B((7=) (2.12) 

For 13 —)■ - 1 - 00 , the above equation yields E{ts{Q)) = 1 — 2E{U'^) (cf. Remark 2.4 of |2]). 

Remark 2.5 A result analogous to that of Theorem 12.21 can be also obtained for x 7 ^ 0; 
however, the involved calculations are very heavy and cumbersome, so we have chosen not 
to develop them. 


Remark 2.6 Once the pair {'g,/3) has been found, such that it verihes fl2.6p . or fl2.7p . it 
may be that 'g is not the Laplace transform of the density function of a random variable U 
with support (O,^). In this case, a solution to the IFPT problem does not exist. This is 
the reason why Theorem 12.21 is formulated in a conditional form. This kind of difficulty in 
showing the existence of a solution to an inverse FPT problem is common to other types 
of inverse problem (see e.g. n. 0,0. a, 0,1221) ; as far as the present IFPT problem 
is concerned, the difficulties are far stronger, because the relation between 'g and / is more 
complicated. For instance, if p = 0 and / is the exponential density with parameter 1, then 
the solution to the IFPT problem with S = 1 and g symmetric in (0,1) does not exist. In 
fact, suppose that the solution exists, then by fl2.7p it follows that 




A2 + 2/3(1 + e-^) 
/3(4 + A2) 


which is not the Laplace transform of a probability density in (0,1), since the third moment 
is negative. 


Taking into account Remark 12.61 we will prove the existence of the density g oi U for a 
class of FPT densities /. For the sake of simplicity, we limit ourselves to the case when 
p = 0, X = 0, S = 1, and g is required to be a function with support in (0,1), which is 
symmetric with respect to the middle point 1/2 ; in fact, for /x 7 ^ 0 the calculations involved 
are far more complicated. 

For an integer k > 0, set /^(A) = f_^e~^'^x^dx; as easily seen, /o(A) = 2sinh(A)/A and the 

recursive relation /^(A) = -h |Jfc_i(A) allows to calculate /^(A), for every k. 

The following Proposition gives a sufficient condition, so that there exists the solution 
{g-,13) to the IFPT problem for the process X with holding and jumps from zero, associated 
to X{t) = Bt, and the barrier S' = 1. 


Proposition 2.7 Let X{t) = Bt, and suppose that the Laplace transform of f{t) has the 
form: 


m = f2kW = 



^,/2/\smh. A/2 

— hk 

(VV 2 ) 


1 + e-^ + A/6 - (1 + ^)e-V^ 

-\/2/Asinh ^a/A/2 

— hk (^\/ A/2j 


(2.13) 

forb > 0 and some integer k > 0. Then, there exists the solution {g,(3) of the IFPT problem 
for X, relative to the barrier S = 1 and the FPT density /, and it results (3 = b and: 

g{u) = g 2 k{u) = (^1 + (1 - {2u - l)^^) , fc > 0, u e (0,1). (2.14) 
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Proof. A simple calculation shows that 



^sinh(A/2) - hki^/^) 


Since g 2 k is symmetric with respect to 5/2 = 1/2, the result follows by inserting g 2 k into 

( 1 ^ . 

□ 

Notice that, letting h go to +oo, fl2.13p becomes equation (2.17) of [2]. 


Remark 2.8 A straightforward calculation shows that, if U E (0,1) has density g 2 k, then 
E{U^) = g(2fc+3) • Then, by using fl2.12p . we obtain that the FPT-distribution corresponding 


to f 2 k has mean E{ti{0)) = | • 


3 Diffusions conjugated to Brownian motion 

In certain cases, a one-dimensional diffusion X can be reduced to BM by a variable change; by 
using this approach, we shall extend to a general one-dimensional DHJ X, which is associated 
to X, the results obtained for Wiener process. For a > 0, let J = [0,a] or J = [0,a), with 
a < -I-CX3, and suppose that X is a time-homogeneous diffusion in J which is the solution of 
the stochastic differential equation (SDE): 

dX{t) = ii{X{t))dt + a{X{t))dBt , X(0) = x E J, (3.1) 

where the coefficients p(x) and a{x) are regular enough functions (see e.g. [2]), so that a 
unique strong solution exists. We consider the following: 

Definition 3.1 We say that X is conjugated to BM if there exists an increasing function 

V : J —;■ IR with n(0) = 0, such that: 

(i) v{x) is continuous for any x E J and it is differentiable in the interior of J; 

(a) v~^{y) is differentiable in the interior of vfJ) and it possesses the right derivative at 
y = 0 and the left derivative at y = v{a), if a < -fcx). 

(Hi) X{t) = v~^{Bt + v{x)), for any t >0. 

Let us suppose that X is conjugated to BM, and consider the associated process X with 
holding and jumps from 0. Notice that holding and random reflection of X at zero corresponds 
to holding and random reflection of Bt + v(x) at zero; moreover the hrst passage ts(x) of 
X through the barrier S, with 0 < 5 < a, corresponds to the hrst passage of Bt -|- v{x) 
through v{S), and ts{x) = Ly5)('i^(a^)), where the superscript B refers to BM. Furthermore, 
let g{u) be the density of the position U E (0, S) from which, once X has hit 0, it starts 
afresh, after the holding time, and let q{y) be the corresponding density of the position 

V = v{U) E (0,n(5)) from which, once Bt + v(x) has hit 0, it starts afresh, after the holding 
time; then, g(y) = g(v~^(y))(v~^)'(y), y E (0,n(5)). Thus, if X is conjugated to BM via the 
function v, then the solution {g,l3) to the IFPT problem for the process X associated to X, 
relative to the FPT density / and the barrier S, can be written in terms of the solution (g, fd) 
to the IFPT problem for the process associated to v{x) -1- Bt, relative to the FPT density 
/ and the barrier v{S), by using that g{x) = q{v{x))v\x). As easily seen, if a; = 0 , g is 
obtained by fl2.6l) with /a = 0, with g in place of g and v{S) in place of S. 
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4 A few examples 


Example 1 (when g is the uniform density in (0, S')) 
Let X{t) = Bt and S' > 0, 6 > 0 and let 


/(A) = 


6(1 - e-^^) 


SVn[\ + 6(1 + e--5^)] - 6(1 - e-^^) 


Then, the solution to the IFPT problem for X relative to S is the pair {g,/3) with (3 = b 
and g{u) = ^1(o,s)(m), i.e. the uniform density in (0, S'); this is easily obtained by searching 
for a solution which is symmetric with respect to S'/2; by using fl2.7p with /S = 6 , we get 
g{\) = ^—§ 3 ;— which is the Laplace transform of g{u) = ^ 1 (o, 5 )(m). In the case S = 1, we can 
obtain the same result by letting k go to inhnity in f 2 k and g 2 k of Proposition 12.71 moreover, 
the mean of the FPT-distribution corresponding to / is 1/3 + 1/6, as it also follows by 

calculating lim^^oo ^ l( 2 k+l) i) Remark 12.81 ). For 6 —)■ +oo (no holding at 0), one 
obtains the function /(A) of Example 3 of P]. 

Example 2 Let X{t) = Bt and S' > 0, 6 > 0 and let 

67r2(l + e-®^) 


/(A) = 


(4AS'^ + 27r2)(A + 6(1 + ’ 


then the solution to the IFPT problem for X relative to S is the pair {g,/3) with [3 = b and 
g{u) = ^ sin (|m) , m e (0, S'). 

In fact, we search for a solution which is symmetric with respect to S'/2; by using fl2.7p 
with /S = 6 , it follows that 'g{Xj = 2 A 2 52+7r2 ’ which is indeed the Laplace transform of 

the function g{u) above. If there is no holding at 0 (i.e. 6 = +cxo), /(A) becomes that of 
Example 1 of [2]. 

Example 3 (when is a Beta density) 

Let X{t) = Bt and S' > 0, 6 > 0 and let 


/(A) = 


6(e-sP^(^^ + 2) + - 2) 


^3A3(1 + e-s^) - 6 (e-^v^(^\/^ + 2) + - 2) + A /6 


Then, the solution to the IFPT problem for X relative to S is the pair {g, /3) with [3 = b and 
g{u) = ^ u{S — u), u E (0, S'). In fact, by using fl2.7p with (3 = b, it follows that 'g{X) = 
^^[e“'^'^(S'A + 2) + S'A — 2], which is indeed the Laplace transform of g{u) = ^ u{S — u). 
Notice that, for S' = 1, is the density g 2 k of Proposition 12.71 for k = 1. For 6 —)■ +cxo (no 
holding at 0), one obtains the function /(A) of Example 4 of [2]. 

Example 4 (when g is the triangular density in (0,1)) 

Let X{t) = Bt and S = 1, 6 > 0 and let 


/(A) = 


26(1 


_ e\/V2^2 


A[A + 6(1 + e“^)] — 26(1 — e 


-Jkj2\2 
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Then, the solution to the IFPT problem for X relative to S is the pair {g,(3) where /3 
and g is the triangular density in (0,1) : 


4m, m e (0,1 
4(1 -m), u e (i,l) 


b 


In fact, by using fl2.7p with /3 = 6, it follows that ^(A) = ^(1 — e V 2 ^ 2 ^ which is indeed the 
Laplace transform of the function g{u) above. For b +oo (no holding at 0), one obtains 
the function /(A) of Example 5 of [2]. 


Example 5. 

A class of diffusions conjugated to BM is given by processes X(t) which are solutions of 
SDEs such as 

dX{t) = -a{X{t))a’{X{t))dt + a{X{t))dBt, X(0) = x (4.1) 

with (t(-) > 0. Indeed, if the integral v{z) = convergent for every z, by Ito’s 

formula, we obtain that X(t) = v~^(Bt + v(x)). Then, by using the results of Section 2, 
examples of solutions to IFPT problems with x = 0 can be easily derived from Examples 1 
to 4, with regard to the process X(t) with holding and jumps from 0, associated to X which 
is driven by the SDE (14.11) . for some choice of a(-); in fact, it suffices to replace S with v(S). 
For instance, diffusions X of this kind are the well-known Feller process (also known as the 
Cox-Ingersoll-Ross (CIR) model), and the Wright & Fisher-like process (see Examples (i) 
and (ii) of | 2 ]). 

Example 6 (Ornstein-Uhlenbeck process) 

Let X (t) be the solution of the SDE: 

dX{t) = —jjiX{t)dt + adBt, X(0) = x, 

where /i,a^are positive constants. By using a time-change, the explicit solution assumes 
the form X{t) = (x -|- B{p{t))) , where p{t) = ^ (e^^* — 1). If S{t) is a moving barrier, 

the FPT of X{t) over S{t) is Tsp) = inf{t > 0 : x -t- B{p{t)) > and so p{js{t)) = 

inf{M > 0 : X -I- > S'(m)}, where S{u) = ^S{p~^{u)). Therefore, if S{t) = 

the IFPT problem for the associated DHJ X(f), and relative to the moving barrier S{t) and 
the FPT distribution E, is reduced to the IFPT problem for DHJ associated to BM, starting 
from X and relative to the constant barrier S'o and the FPT distribution F = F o p~^. Thus, 
if X = 0, explicit examples for the Ornstein-Uhlenbeck process and the exponential barrier 
S(t) = S'oe“^*, can be easily derived from Examples 1 to 4. 

Example 7 (Geometric Brownian motion) 

Let X (t) be the solution of the SDE: 

dX{t) = rX{t)dt -|- aX{t)dBt, X(0) = x > 0, 

where r and a are positive constant. This is a well-known equation in the framework of 
Mathematical Finance, since it describes the time evolution of a stock price X. The explicit 
solution is X{t) = where p = r — cr^/2. Let us consider the moving barrier S{t) = 

gO-s'+Mt. the IFPT problem for the associated DHJ X(t), relative to S{t) and the 

FPT distribution E, is reduced to the IFPT problem for DHJ associated to BM with drift 
(/i — p')/(J, starting from ^ and relative to the constant boundary S, and the same FPT 
distribution E. 
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